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Abstract 

We present the applications of variational-wavelet ap- 
proach to nonlinear (rational) model for spin-orbital mo- 
tion: orbital dynamics and Thomas-BMT equations for 
classical spin vector. We represent the solution of this dy- 
namical system in framework of periodical wavelets via 
variational approach and multiresolution. 

1 INTRODUCTION 

In this paper we consider the applications of a new numeri- 
cal-analytical technique which is based on the methods of 
local nonlinear harmonic analysis or wavelet analysis to the 
spin orbital motion. Wavelet analysis is a relatively novel 
set of mathematical methods, which gives us a possibil- 
ity to work with well-localized bases in functional spaces 
and give for the general type of operators (differential, inte- 
gral, pseudodifferential) in such bases the maximum sparse 
forms. Our approach in this paper is based on the general- 
ization of variational-wavelet approach from [l]-[8], which 
allows us to consider not only polynomial but rational type 
of nonlinearities [9]. The solution has the following form 
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which corresponds to the full multiresolution expansion in 
all time scales. Formula (^ gives us expansion into a slow 
part z^""^ and fast oscillating parts for arbitrary N. So, we 
may move from coarse scales of resolution to the finest one 
for obtaining more detailed information about our dynami- 
cal process. The first term in the RHS of equation (1) corre- 
sponds on the global level of function space decomposition 
to resolution space and the second one to detail space. In 
this way we give contribution to our full solution from each 
scale of resolution or each time scale. The same is correct 
for the contribution to power spectral density (energy spec- 
trum): we can take into account contributions from each 
level/scale of resolution. 

In part 2 we consider spin-orbital motion. In part 3 start- 
ing from variational formulation we construct via multires- 
olution analysis explicit representation for all dynamical 
variables in the base of compactly supported periodized 
wavelets. In part 4 we consider results of numerical cal- 
culations. 

2 SPIN-ORBITAL MOTION 

Let us consider the system of equations for orbital motion 
and Thomas-BMT equation for classical spin vector [10]: 
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dq/dt = dHorb/dp, dp/dt = —dHorb/dq, ds/dt 
w X s, where 
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q = {qi,q2,q3),p = (pi,P2,P3) are canonical position 
and momentum, s — (si, ,S2, S3) is the classical spin vec- 
tor of length fi/2, tt = (tti, 7r2, tts) is kinetic momen- 
tum vector We may introduce in 9-dimensional phase 
space z ~ {q,p,s) the Poisson brackets {f{z),g{z)} ~ 
fq9p ~ fp9q + [fs X 5s] ■ s and the Hamiltonian equations 
are dz/dt = {z, H} with Hamiltonian 



H = Horbiq,P,t) + w{q,p,t) ■ s. 



(3) 



More explicitly we have 
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We will consider this dynamical system in [11] via invari- 
ant approach, based on consideration of Lie-Poison struc- 
tures on semidirect products. But from the point of view 
which we used in [9] we may consider the similar approx- 
imations and then we also arrive to some type of polyno- 
mial/rational dynamics. 

3 VARIATIONAL WAVELET APPROACH 
FOR PERIODIC TRAJECTORIES 

We start with extension of our approach to the case of peri- 
odic trajectories. The equations of motion corresponding 
to our problems may be formulated as a particular case 
of the general system of ordinary differential equations 
dxi/dt — fi{xj,t), j = l,...,7i), < t < 1, where 
fi are not more than rational functions of dynamical vari- 
ables Xj and have arbitrary dependence of time but with 
periodic boundary conditions. According to our variational 
approach we have the solution in the following form 

=a;,(0) + ^Af(^fe(i), x,(0) = x,(l), (5) 



where are the roots of reduced algebraical systems of 
equations with the same degree of nonlinearity and ipk (t) 
corresponds to useful type of wavelet bases (frames). It 
should be noted that coefficients of reduced algebraical sys- 
tem are the solutions of additional Unear problem and also 
depend on particular type of wavelet construction and type 
of bases. 

Our constructions are based on multiresolution appro- 
ach. Because affine group of translation and dilations is 
inside the approach, this method resembles the action of 
a microscope. We have contribution to final result from 
each scale of resolution from the whole infinite scale of 
spaces. More exactly, the closed subspace Vj {j G Z) cor- 
responds to level j of resolution, or to scale j. We consider 
a r-regular multiresolution analysis of _L^(R") (of course, 
we may consider any different functional space) which is a 
sequence of increasing closed subspaces Vj: 

...V-2 C V-i C C C F2 C ... (6) 

Then just as Vj is spanned by dilation and translations of 
the scaling function, so Wj are spanned by translations and 
dilation of the mother wavelet tpjk (x), where 

^Pjkix)=2i/^^p{2^x-k). (7) 

All expansions, which we used, are based on the following 
properties: 
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L^{R) = Vo^Wj (8) 
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We need also to find in general situation objects 

00 

^tt::t= J n'^*(^)d^' (9) 

but now in the case of periodic boundary conditions. Now 
we consider the procedure of their calculations in case 
of periodic boundary conditions in the base of periodic 
wavelet functions on the interval [0,1] and corresponding 
expansion (1) inside our variational approach. Periodiza- 
tion procedure gives us 

'fjA^) = ^^j,k{^ - ^) (10) 

So, (p, tj) are periodic functions on the interval [0,1]. Be- 
cause (fij^k = <fij,k' ii k = fc'mod(2-'), we may consider 
only < k < 2^ and as consequence our multiresolution 
has the form [J V^" = L^[0, 1] with Vj = span{(^j,/j;}^-ij^^ 

3>0 

[12]. Integration by parts and periodicity gives useful 
relations between objects (9) in particular quadratic case 
{d = di+ d2): 
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So, any 2- tuple can be represented by A^. Then our sec- 
ond additional linear problem is reduced to the eigenvalue 
problem for {A^}o</c<2j by creating a system of 2^ ho- 
mogeneous relations in A^ and inhomogeneous equations. 
So, if we have dilation equation in the form ip{x) = 
\/2 X^feez hk^p{2x — k), then we have the following ho- 
mogeneous relations 

^fe = 2' E E >^mhA^,^,,_^, (12) 
m=0 e=o 

or in such form AA<* = 2''A'^, where A"^ = {Ai}o<k<2i- 
Inhomogeneous equations are: 

Y^Mfh.i = d\2-i''^, (13) 

I 

where objects M/(|£| < N — 2) can be computed by re- 
cursive procedure 

= 2-^'(2'^+i)/2^//, (14) 
Mfe =< x\ ^o,e >= E C") "'^"'^0. ^0 = 1- 

So, we reduced our last problem to standard linear alge- 
braical problem. Then we use the methods from [9]. As a 
result we obtained for closed trajectories of orbital dynam- 
ics the explicit time solution (1) in the base of periodized 
wavelets (10). 
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Figure 1: Periodic wavelet 
4 NUIVIERICAL CALCULATIONS 

In this part we consider numerical illustrations of previous 
analytical approach. Our numerical calculations are based 
on periodic compactly supported Daubechies wavelets and 
related wavelet families (Fig. 1). Also in our modelling we 
added noise as perturbation to our spin orbit configurations. 

On Fig. 2 we present according to formulae (2), (6) con- 
tributions to approximation of our dynamical evolution (top 
row on the Fig. 3) starting from the coarse approximation, 
corresponding to scale 2° (bottom row) to the finest one 
corresponding to the scales from 2^ to 2^ or from slow 



to fast components (5 frequencies) as details for approxi- 
mation. Then on Fig. 3, from bottom to top, we demon- 
strate the summation of contributions from corresponding 
levels of resolution given on Fig. 2 and as result we re- 
store via 5 scales (frequencies) approximation our dynam- 
ical process(top row on Fig. 3 ). The same decomposi- 
tion/approximation we produce also on the level of power 
spectral density in the process with noise (Fig. 4). 
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Figure 2: Contributions to approximation: from scale 2^ to 
2^ (with noise). 
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Figure 3: Approximations: from scale 2^ to 2^ (with 
noise). 
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Figure 4: Power spectral density: from scale 2^ to 2^ (with 
noise) 
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